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Abstract 

The phase transitions between the XY, the dimer and the Haldane 
phases of the spin-1 bond-alternating XXZ chain are studied by the nu- 
merical diagonalization. We determine the phase diagram at T = and 
also identify the universality class with the level spectroscopy. We find 
that exactly on A = there is a Berezinskii-Kosterlitz-Thouless transi- 
tion line which separates the XY and the Haldane phase and there exists a 
multicritical point of the XY, the dimer and the Haldane phases. We dis- 
cuss that the critical properties of this model is of the 2-D Ashkin- Teller 
type reflecting the hidden Z2 x Z2 symmetry. 

Haldane predicted the difference of the critical properties between integer 
S and half odd integer S spin chains. When S is an integer, the Heisenberg spin 
chain has a unique disordered ground state with the finite energy gap, while 
when S" is a half odd integer it has no energy gap and belongs to the same 
universality class as the S = 1/2 case. Den Nijs and Rommelse |^ argued that 
the Haldane gap system is characterized by the string order parameter which is 
defined by the non-local operator. Recently the symmetry breaking is recognized 
as the breakdown of the hidden Z2 x Z2 symmetry In this letter we study 
the phase diagram of the 5 = 1 anisotropic spin chain with bond-alternation 

^ = E(i + s{-iy)[s^s^+, + + A5|5|+i], (1) 

to see this symmetry breaking. By studying the phase diagram, we show the 
difference of the topology of the phase diagrams between 5 = 1 and 5=1/2 
cases. 

First we review the two dimensional (2-D) Ashkin- Teller model which will 
be useful to understand the model (^. The 2-D Ashkin- Teller model can be 
thought as two Ising models coupled by a four-spin interaction, and has the 
Z2 X Z2 symmetry. It is known that there is the 2-D Gaussian critical line of 
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continuously varying critical exponents Q|, and at one end, this line breaks up 
into two 2-D Ising critical lines At another end of the Gaussian critical 
line, it meets two Berezinskii-Kosterlitz-Thouless (BKT) critical lines and 
so there is a massless region (so called "critical fan") in the 2-D Ashkin- Teller 
model. The 2-D Gaussian critical line separates the fully ordered phase and 
the fully disordered phase about the Z-i x Z-i symmetry, and the two 2-D Ising 
critical lines are the boundaries of the partially ordered phase. According to 
Kohmoto et al. the 2-D Ashkin- Teller model can be mapped to the S" = 1/2 
XXZ chain with bond- alternation whose Hamiltonian is defined by eq.(|^). For 
this S' = 1/2 model the Gaussian critical line corresponds to the (5 = line 
(— 1/\/2 < A < 1) and at the A = 1, 5 = point the bifurcation to the 2- 
D Ising critical lines occurs. In finite systems, it is known that the boundary 
condition is important for the eigenvalue structure of these models 0. For 
the periodic boundary conditions, the 5* = 1/2 XXZ chain is ?7(1) symmetric 
and the quantum Ashkin- Teller model is Zi x and the correspondence of 
the eigenvalues for both models is not complete. On the other hand, for open 
boundary conditions the eigenvalues of the = 1/2 XXZ chain with 2M sites 
are exactly related to those of the self-dual M-site quantum Ashkin- Teller chain. 

Let us return to the S = \ case. For (5 = case this model has been studied 
numerically |l^ |ll|, |l^, ^ in relation to the Haldane conjecture In the 
region A > Ac2 = 1.16 ± 0.02 the system is Neel ordered and it has a two-fold 
degenerate ground state. When Ac2 > A > Ad ~ 0.0, the system is in the 
Haldane phase, where the ground state is a unique singlet with an energy gap 
and the spin correlation decays exponentially. In the region — 1 < A < Ad 
the system is the XY phase, characterized by the power law decay of the spin 
correlation hmction and the gapless excitation spectrum. For A < — 1 this 
system is ferromagnetic long-range ordered. The phase transition at Ac2 belongs 
to the 2-D Ising universality class ||ll|, ^2) . It is thought that the BKT transition 
occurs at Ad. Numerically Botet and JuUian estimated Ad '--^ 0.1, Kubo and 
Takada Ad ~ 0.3, Sakai and Takahashi Ad = —0.01 ± 0.03 and Yajima and 
Takahashi Ad = 0.069 ± 0.003. However it is difhcult to evaluate the precise 
value of the Ad from the gap behavior of finite size systems because of the 
singular behavior of the excitation gap ^ exp(— const/\/A — Ad). In fact, 
the finite-size scaling method may lead to false conclusions for the BKT-type 
transition |l^, So in this letter, we mainly discuss the method to obtain 
the BKT transition points. 

From the other point of view, the isotropic case with bond-alternation 

H = ^(l + (5(-l)^)S,.5,+i (2) 

3 

has been extensively studied. Theoretically Affleck and Haldane ||l^, |l^ studied 
the Heisenberg case for the arbitrary S*, mapping to the 0(3) nonlinear cr-model, 
and showed that the topological angle is given by 6' = 27r5(l + S) and the 
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system should be massless when Q jl-K is half odd integer. Numerically several 
authors estimated the transition point be between the Haldane phase and the 
dimer phase, and calculated the conformal anomaly c using the density matrix 
renormalization group method |2^] , the Monte Carlo method |^ and the exact 
diagonalization |2^. They concluded 6^ — 0.25 and c = 1. Recently Totsuka 
et al. calculated the critical dimensions up to the logarithmic corrections 
and identified the universality class as the k — 1 SU{2) Wess-Zumino-Witten 
(WZW) model. 

According to Kennedy and Tasaki in the Haldane phase the hidden 
Z2 X Z2 symmetry is fully broken, in the Neel phase only one Z2 symmetry is 
broken and in the dimer phase there is no symmetry breaking. The edge states 
with open boundary conditions reflect this symmetry. This suggests that the 
transition point at A = 1, (5 ~ 0.25 is the Ashkin- Teller bifurcation point. 

On the anisotropic case with bond alternation, Singh and Gelfand 
treated this model by series expansions. They stated that there exists a mul- 
ticritical point of the Neel, the dimer, the Haldane and the XY phases. In 
their phase diagram there is no 2-D Gaussian critical line, so it contradicts to 
the expectation that the transition point A = 1, ^ ~ 0.25 is the Ashkin- Teller 
bifurcation point. 

In the following we consider the XY-dimer and the XY-Haldane phase tran- 
sitions. Using the bosonization method, Schulz showed that the effective 
Hamiltonian of this transition is described by the sine-Gordon model: 



H= — I dx 
2tt 



2 , V fdcj)'^'^ 
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dx cos \/8(/i, (3) 



(27ra)2 

where H is the momentum density conjugate to (p, 

[cj){x),n{x)]^i5{x~x), (4) 

a is the lattice constant and v is the spin wave velocity. The dual field 6{x) is 
defined as dx6{x) = 7rH(a:). Here we neglect several irrelevant terms for the BKT 
transition. In the full continuum Hamiltonian there is another part which is a 
continuum limit of the 2-D Ising model and it is important for the Neel-dimer 
transition and the Neel-Haldane transition. The BKT transition occurs when 
the scaling dimension of the operator cos y/Scj) becomes 2 or the rcnormalizcd K 
is 1. 



According to Cardy 1 25 , assuming conformal invariance of correlation func- 
tions, the relation between the critical dimension and the energy gap of the 
finite system with periodic boundary conditions is 

E,,{L)-Eo{L)^^^ (5) 

where v is the spin wave velocity of the model, Xn is the scaling dimension and L 
is the number of the unit cell (for the present model L = N/2, N is the number 
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of sites). We determine the spin wave velocity from the energy gap at the wave 
number k = 2tt/L and J2 = ^ as 

LAE(k = 27r/L) 
V = hm v{L) = lim ^- (6) 

Numerically we fitted values of v{L) to the form v{L) = i; + a/L + b/L'^. The 
leading finite size correction of the ground state energy is 

Eo{L) = eoL-—, (7) 

where c is the conformal anomaly number. 

Besides the operator cos -s/St/), there are several operators whose scaling di- 
mensions become 2 at the BKT transition point. Recently one of the authors 
[ p4| studied the structure of these operators for finite size systems. He found 
that at the BKT transition point cos ^/8(f) and the marginal operator 

which is identical to the free part of the Lagrangian density, should be hybridized 
to satisfy the orthogonality condition of the two-point correlation functions in 
the fields. This hybridization comes from the existence of cosVScj) term in 
eq.(^) which generates the wave function renormalization. He showed that the 
ratio of leading logarithmic corrections of the scaling dimensions of cos y/Scjylike, 
sm\/8(f>, exp{^iy/89) and marginal-like operators (in the following we define 
the scaling dimensions of these operators as xi, X2, X3 and xq respectively) 
is 2 : 1 : — 1 : — 1. Near the transition point the scaling dimensions xq and 
X3 cross linearly. Note that in the XY phase the scaling dimension of the 
marginal operator M remains 2 but for other operators the scaling dimensions 
vary continuously as the effective K varies. Since the scaling dimension is related 
to the excitation energy as eq.(^, this means that at the critical point there 
exists a degeneracy of excited states which correspond to the exp{^iy/89) and 
the marginal-like operator. So we can determine the BKT transition point by 
the level crossing of these excitations, and we can also eliminate the logarithmic 
corrections by the appropriate average of the excitations. We call this method 
the level spectroscopy. 

In order to identify the excitations with those of the sine-Gordon model, 
we use the following symmetries. Hamiltonian (|l|) is invariant under spin ro- 
tation around the z-axis, translation by two-site (Sf —^ S'f^.j) , space inversion 
{Sf ^N-i+i) ^"^^ ^P^^ reversal (S'f —Sf). Corresponding eigenvalues are 
J2iLi q — 4:TTn/N, P — ±1 and T — ±1 respectively. The symmetry opera- 
tions in the sine-Gordon model are as follows |24j. The operation to the space 
inversion (P) is 

0^-0, 9^9 + V2Tr, X ~* -X, (9) 
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and the operation to the spin reverse (T) is 

4>^-4>, e^-e + V2Tr. (lo) 

The eigenvalues of these symmetry operations are summarized in Table 1. For 
finite systems the ground state is singlet = 0, q = 0, P = T = 1) in the 

whole region (A > —1). In ^ S*^ = 0, g = subspace, the first excited state is 
P = T ~ 1, the second one is P = T = —1 and the third one is P = T = 1 in the 
masslcss phase. In ^ S'^ = ±4, g = subspace the first excited state is P = 1. 
On the BKT line the first ^ S*^ = 0, g = 0, P = T = 1 and -5"^ = ±4, g = 0, 
P = 1 excitations degenerate. Note that in ^ 5^ = ±1, P = —1 subspace the 
lowest excitation states correspond to the operators exp(=Fz6'/\/2) in the sine- 
Gordon language, whose scaling dimension is x — 1/8 {rj — 1/4) at the BKT 
transition point. 

In numerical calculations, we consider finite rings of N sites {N = 6, 8, 10, 12) 
with periodic boundary conditions. We calculate the excitation energies of 
X] 5*^ = and ^ S*^ = ±4 corresponding to cos VStj), sin VStj), marginal opera- 
tor and exp(=Fi-\/80). In Fig. 1(a) we show the excitation energies for iV = 12, 
A = —0.5. It is seen from this figure that the marginal-like excitation = 
0) and the excitation corresponding to the operator exp(=Fi-\/80) = ±4) 

cross linearly. We determine the critical point of the XY-dimer phase transi- 
tion by the crossing point of these excitations as mentioned above. Figure 1(b) 
shows the excitation energies for = 12, A = 0. In this case one of the state 
(^ 5-^ = 0, g = 0, P = T = 1) and the eigenstate (^ S*^ = ±4, g = 0, P = 1) 
are degenerate in the whole region (numerically this degeneracy is exact even for 
finite systems), thus we conclude that there is a BKT transition line between the 
XY and the Haldane phases on A = 0. Finally we obtain the multicritical point 
of the XY, the dimer and the Haldane phases by the crossing of the excitations 
corresponding to the operators sini/St/) — 0,P = —1) and exp(=Fi'y86') 

{J2 = ±4, P = 1). This multicritical point corresponds io K — 1 and g — Q 
in the continuum Hamiltonian (^). We estimate this point as 5 — 0.22 ± 0.01. 
Figure 2 indicates the full phase diagram of this model. Other boundaries (be- 
tween the dimer, the Haldane and the Neel phases) are determined by the usual 
finite-size scaling method. 

In order to verify that the XY-dimer and the XY-Haldane phase transitions 
are of the BKT-type, we present (xq -I- a;2)/2 and (2x0 + xi)/2> in Fig. 3 which 
eliminate the logarithmic corrections of the scaling dimension. Since we only 
investigate small size systems, these values are not exactly equivalent to 2, due 
to the existence of irrelevant fields, but they converge rapidly to this value 
with increasing system size. Also we estimate the conformal anomaly using 
eq.(0), as c = 0.97 for A = -0.5 and c = 0.98 for A = 0, (5 = 0. Hence we can 
conclude that the transitions are of the BKT-type. On the other hand the phase 
transition between the Haldane and the dimer phases is thought to belong to 
the 2-D Gaussian universality class. In this case the operator cos in the 
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effective Hamiltonian eq.(|^) is relevant, that is, the effective K is smaller than 
1 and the second order transition occurs at g 

Now we consider the reason why the boundary of the XY and the Haldane 
phases is just on the A = Hne. On this Hne ^ 5*^ = and ^ 5^ = ±4 
excitations are degenerate exactly. Note that in the case of the open boundary 
conditions, similar degeneracy was reported by Alcaraz and Moreo |^6|, who 
discussed that it strongly suggests A^i = 0. Although we have no explicit 
explanation on these degeneracies, we remark the following symmetry. By the 
unitary operator U = exp[j7r^^-^(j + \)S^\, the Hamiltonian is transformed 
as UH{A)U^^ ~ —H{—A). This means that there is an additional symmetry 
at A = 0. When A = 0, for an eigenstate \tpo), U\^po) is also an eigenstate 
of H. Perhaps this symmetry makes the phase boundary of the XY and the 
Haldane phases on the line A = for the S = 1 case. Relating to this, Alcaraz 
and Hatsugai |^ showed that the S* = 1 XY model does not have the string 
order, and their argument can also be applied to the bond-alternation cases. In 
comparison, the S = 1/2 XY model with bond-alternation is equivalent to the 
free fermion model. In this case, there exists the degeneracy between ^ S*^ = 
and ^ 5^ = ±2, and the second order transition occurs at (5 = 0. 

Let us discuss the relation between the present phase diagram and the 
bosonization by Schulz. The renormalization flow of eq.(^) is shown in Fig. 4. 
Our phase diagram Fig. 2 means that all the starting points of the renormaliza- 
tion for A = and < (5 < 5c he just on the line in Fig. 4 which separates the 
XY and the Haldane phases, and the multicritical point of the XY, dimer, Hal- 
dane phases ( A = 0, 5 = Sc) corresponds to the origin of Fig. 4. Although Schulz 
stated that the bare values of K and g are K = a/I -I- (2A/7r) (independent of 
6) and g = — A + (5^ respectively, it is easy to see that his expression cannot ex- 
plain the above-mentioned behavior. We consider that this discrepancy comes 
from that Schulz dropped many operators in the course of the bosonization. 
However, apart from the bosonization, the BKT multicritical structure can be 
expected based on the U{1) x Z2 x Z2 symmetry [Q. If the bare values have 
the form, for instance, K = -^/l -I- (2A/7r) — aiS'^ and g = 02 — A -t- a^S^ with 
positive constants ai, 02 and 03, the renormalization behavior of eq.(||) is con- 
sistent with our phase diagram Fig. 2. Unfortunately we could not derive such 
expressions from the original spin Hamiltonian. To derive the correct form of 
the effective Hamiltonian is a future problem. 

Lastly, we argue the difference of the phase diagram between the 5 = 1/2 
and the S = 1 cases. Note that the phase diagram is symmetric on the trans- 
formation S —(5. Reflecting the hidden Z2 x Z2 symmetry, the S" = 1/2 
model and the S = 1 model with bond alternation have the same phase di- 
agram with the 2-D Ashkin- Teller model. For the S' = 1/2 model, there are 
two BKT lines, one 2-D Gaussian critical line and two 2-D Ising type critical 
lines for —1 < i5 < 1. On the other hand, the = 1 model has three BKT 
lines, two 2-D Gaussian lines and three 2-D Ising type critical lines. Oshikawa 
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[ p8| studied the hidden Z2 x Zi symmetry of the arbitrary integer S Heisenberg 
Hamiltonian with bond-alternation (|^). In the content of the VBS picture, we 
expect the successive dimerization transition which is in agreement with the 
Afheck-Haldane prediction. From these considerations, we can expect that for 
the arbitrary S XXZ chain with bond-alternation, there are 25-1-1 BKT lines, 
25* 2-D Gaussian lines and 25 4- 1 2-D Ising lines. In addition, for half odd 
integer 5, there is a Gaussian line in the region (5 = 0, — 1 < A < 1 and at 
A = 1 the model becomes the fc = 1 5C/(2) WZW model. 

In summary, we present the phase diagram of the 5=1 bond-alternating 
XXZ chain (Fig. 2) which is determined by the level crossing of some special 
excitations on the basis of renormalization group analysis (Fig.l) and also show 
the logarithmic corrections of the critical dimensions on the BKT critical line 
(Fig. 3). There exists a BKT line which separates the XY and the Haldane 
phases on A = line where the X] "S"^ = and the X] '5'^ = ±4 eigenstates 
are degenerate on the whole line. The critical properties are recognized as that 
of the 2-D Ashkin- Teller model. More detailed studies of this model will be 
presented in elsewhere. 
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Symmetry operation 
q T P 


Operators of 
sine-Gordon model 


Notation of 
scaling dimensions 





1 1 


M 


.To 





1 1 


cos ^/S4' 


Xi 





-1 -1 


sin\/8(^ 


Xl 


±4 


* 1 


exp =Fi\/8^ 





Table 1: Eigenvalues of symmetry operations 
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Figure Caption 

Fig.l Excitation energies with q = OX -S*^ = and 5''' ±4 for (a) N 12, 
A = -0.5 and (b) = 12, A = 0. O's arc ^ S*^ = 0, g = 0, P = 1, T = 1 
excitations, □ is ^ S'^ = 0, g = 0, P = —1, T = —1 excitation and + is the 
^ 5^ = ±4, q = 0, F = 1 excitation. 

Fig. 2 The phase diagram in the AS plane. F means the ferromagnetic long- 
range ordered region. The XY-dimer and the XY-Haldane phase boundaries are 
of the BKT type, the dimer-Haldane boundary is of the 2-D Gaussian type and 
the Neel phase boundaries are of the 2-D Ising type. 

Fig. 3 Scaling dimensions eliminating logarithmic corrections (xi + X3)/2 (+) 
and (x2 -I- 2x3) /3 (□) at the BKT transiton points, (a) plots on the BKT line 
between the XY and the dimer phases and (b) plots on the BKT line between 
the XY and the Haldane phases (A = 0). Bare values are also presented (O). 
Fig. 4 The renormalization flow of eq.(|3|). 
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